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Abstract 

We calculate the radiative decay amplitude for Majorana neutrinos in trilinear 
i?-parity violating supersymmetric framework. Our results make no assumption re- 
garding the masses and mixings of fermions and sfermions. The results obtained 
are exemplary for generic models with loop-generated neutrino masses. Compar- 
ison of this amplitude with the neutrino mass matrix shows that the two provide 
independent probes of CP-violating phases. 



1 Introduction 

In the coventional minimal supersymmetric standard model (MSSM) without right- 
handed neutrinos and with conserved i?-parity it is not possible to write renormalizable 
operators which generate neutrino masses. However, conservation of lepton and baryon 
numbers leading to exact i?-parity is an ad-hoc symmetry in the MSSM, and there is 
no deep underlying principle as to why such a discrete symmetry will remain conserved. 
Once we admit i?-parity violating (RPV) couplings P in the superpotential, neutrino 
masses are generated through diagrams involving RPV couplings. These are Majorana 
masses which result as a consequence of broken lepton number (s). There are two kinds of 
explicit RPV couplings: bilinear and trilinear. The bilinear ones induce sneutrino vacuum 
expectation values, allowing neutrinos to mix with the neutralinos, and in this mecha- 
nism only one physical neutrino obtains mass at the tree level |2] . Trilinear Yukawa type 
couplings generate a complete neutrino mass matrix through one-loop self-energy graphs 
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Obviously, one can take a combination of bilinear and trilinear couplings, along with 
possible soft terms, for a complete analysis (4j. Here we concentrate only on the trilinear 
superpotential couplings. Importantly, the diagrams that generate neutrino masses also 
generate dipole moments by insertion of photons to the internal lines. Clearly, in the 
absence of Dirac-type operators, such dipole moments are transitional. 

The connection between neutrino Majorana masses and transition moments have been 
studied in the past jSUHl and predictions have been made using the constraints on trilinear 
RPV couplings placed from the data on neutrino masses and mixings [7 (compare the 
present update in jH]). However, as regards the analytic formulation of dipole moments, 
two crucial things were still missing. One of these is a general analysis with the inclusion 
of fermion and sfermion flavor mixings. The other is a careful handling of the phases from 
different sources that enter the analysis. We try to address these two issues in the present 
paper. In doing so, we ask the following question: how different is the combination of 
phases that appears in the neutrino radiative decay from the one that appears in the 
neutrinoless double beta decay? We keep the formalism at a very general level. Thus 
our results can be applied, with appropriate changes in notation, to a general class of 
models P involving Majorana neutrinos, although for the sake of definiteness we choose 
the supersymmetric RPV model, for which the lepton number violating part of the RPV 
superpotential can be written as 

^ = ^ E >^^^'^L,L,E^, + K,kUQ,Dl . (1) 

ijk ijk 

Here Lj denotes the doublet leptonic superfield and Qi is the doublet quark superfield, 
whereas and denote SU(2)-singlet superfields which contain the left-chiral compo- 
nents of charged antileptons and down-type antiquarks. Stringent upper limits exist on 
all these couplings from different experiments [TIH lllj. 

2 Feynman rules with fermion and sfermion mixings 

Without any loss of generality, we can take the quark superfields in a basis such that the 
up-type quark mass matrix will come out to be diagonal. Also, the lepton superfields 
will be defined in a basis where the charged lepton mass matrices will come out to be 
diagonal. In addition, we can take the singlet fields in a way that they will not have to 
be diagonalized further. This will be called the flavor basis. 

The mass matrices for neutrinos and the d-type quarks will be non-diagonal in this 
basis. We define the mass eigenstates as follows: 

^iL = UiafaL i 

a 

diL = ^ViadaL, (2) 

a 

where in each case, the states on the left side are the flavor states, and the states on the 
right are mass eigenstates. The matrices U and V are respectively the PMNS and the 
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CKM mixing matrices for leptons and quarks. The squark fields will also mix in general. 
We will denote the superpartners of the left-chiral quark fields cLl by without the 
subscript L since the scalar field does not carry any chirality. Similarly d^, the conjugate 
of right-chiral quark fields cLr, have superpartners which will be denoted by d'^. In the 

mass matrix, the fields d will mix with the fields d^ . For three generations, there will be 
then six eigenstates for down squarks. These will be denoted by Ayi, where 

dk = ^ KkA^A , 
A 

dfj- = J^KkA^A. (3) 
A 

Here, for n fermion generations, both K and K are nx2n matrices, which are respectively 
the first and last n rows of a 2n x 2n unitary mixing matrix K of the down squark sector. 

From the superpotential given in Eq. ([Q), when we derive the Lagrangian containing 
ordinary fields, any cubic term in the superpotential produces various kinds of Yukawa 
couplings. We are interested in processes with external neutrino lines, so we seek for terms 
where the neutrino field (and not the sneutrino field) appears. For the time being, let us 
deal with only the second term, i.e., the A' term. We shall comment on the inclusion of 
the A term later. The relevant Yukawa couplings are then: 

-^Y = Y.^^3k^~^LC-'d,J% + Y,K,k^lC-'dlJ,+\i.c.. (4) 

ijk ijk 

Using Eq. this can be rewritten in terms of mass eigenstates in the form 

= Y.^U^IC-'P,daA\ + Y,KaA^IC-'P,d:AA + h.C., (5) 



aaA aaA 



where 



^aaA — ^jk^iaVjaKf^A i 

ijk 

^'aaA = "^KjkUiahaKjA- (6) 

ijk 

The appearance of the Kronecker delta in the latter formula is a reminder that for singlet 
fields, the flavor indices and the mass indices are interchangeable. In Eq. we have 
also used the notation for chiral projection operators that will be used henceforth: 

Pl = ^(1-75), PR = i(l + 75). (7) 

Contributions to radiative neutrino decay occur at the 1-loop level through diagrams 
such as those shown in Fig.^ If the matrix K is block-diagonal, i.e., there is no mixing be- 
tween the fields d and the fields d'^ in the mass matrix, contributions will be proportional 



3 




Figure 1: 1-loop diagrams for the neutrino photon vertex. All fields shown are mass 
eigenstates. The photon attaches to the quark in diagram (a), and to the squark in (b). 

to the external neutrino masses and should therefore be small. We will be interested 
in the more general case where K is not necessarily block-diagonal. Diagrams for this 
general case are shown in Fig. ^ As we will see, in this case the radiative neutrino decay 
amplitude will be proportional to the masses of the down-type quarks. 

Before calculating the amplitude, let us put Eq. (0) is a more conventional form. We 
have used the symbol d!^. This is the conjugate of the (i-quark field. For any fermion field 
if), the Lorentz-covariant conjugate can be defined as 

r = locr ■ (8) 

It is then easy to see that 



Antisymmetry of the matrix C follows from the relation {'il)'^y = ij^. Using this as well as 
the anticommutation of fermion fields, we can write 

^IC-'P^d: = dl^C-'P^U^ = daPLl^a , (10) 

using Eq. in the last step. 

Further, for Majorana neutrino fields, the conjugate field is the same as the original 
fields, apart from possibly a phase factor: 

= rjaU^ . (11) 
Thus, for Majorana neutrino fields, we obtain 

uJC-^ =l^a = Va^a ■ (12) 

Putting Eqs. (fTUI) and (fT^ back into Eq. (jSj), we obtain 

= E [^^{vaKaAPL + A1^Pr)4AJ, + da^vXaA^R + AUPLj^^.A^] , (13) 

aaA 
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Figure 2: Same as in Fig. Q except the internal lines are now conjugated. 



where we have written the hermitian conjugate terms explicitly. By using Eqs. 
(fTTj). these interactions can also be written as 



and 



'^y = Y. [< (^-aPl + <a1aPr) ^.a^ + (aIaPr + ^«AUPl) d^AA] . (14) 



aaA 



This form will be useful when we calculate the diagrams shown in Fig. |21 where the 
internal particles are conjugated. As is well-known, such cojugated diagrams exist when 
we deal with Majorana neutrinos jT2]- 



3 Radiative neutrino decay amplitude 
3.1 Calculation of the vertex function 

The vertex function F'^ will be defined in such a way that the matrix element of the photon 
vertex with on-shell fermions is given by 

M = Mp')T\^{p)Axiq) . (15) 

For the diagram in Fig. the contribution to the vertex function will be denoted 
by r^^, which can be easily written down, using the vertex Feynman rules from the 

interaction given in Eq. (jlHj) . The result is 

,T^^ _ f d^k ( ^, p , p y(^-^ + m,J 

'^B) - 1^ J (2^^1^'^VaPl + A,.^PRj _ _ 



A 



where is the charge of the d-quark. From a general analysis of the form factors related 
to the electromagnetic vertices, we know that the radiative transition amplitudes can 
have only the magnetic and electric dipole form factors. These are helicity flipping terms. 
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containing an even number of Dirac matrices sandwiched between the spinors. We can 
separate out such terms and write 



led 



X 



aA 



{2n) 



(17) 



Similarly, we can write down the contribution of Fig. ^p. By combining the factors in 
the denominator of the integrand, we can identify the magnetic and electric dipole terms, 
which are of the form 



m 



;i8) 



aA 



Contributions to JaA from Fig. [T^ and Fig. [Dd are given by 



J. 



IB) 



aA 



J. 



IB) 



aA 



dx 



dx 



d'^k 



-2(1 - x) 



{2^Y 
d^k 



F - xMl 



3 ' 



-2x{l - x) 



(2^) 



1 3 • 



(19) 
(20) 



k"^ — XM\^ ~" (1 ~ ^)^da 

In writing these integrals, we have neglected neutrino masses in the denominators. Thus 

d^k -2(1 - x) 



J. 



m 

aA 



dx 



i2n) 



21 



F - xMl^ - (1 - x)ml^ 



We now come to the conjugated diagrams of Fig. |2l which are related to the diagrams 
of Fig. [U through 



dn 



d': 



A, 



(22) 



It is now convenient to use the Yukawa couplings in the form given in Eq. (|14j). Com- 
paring this with Eq. (jl3j) . we see that the corresponding amplitudes are related by the 
substitutions 



-^aaA * ^ ^aaA 



(23) 



and an overall negative sign because the photon line now attaches to an internal line of 
opposite charge. So we obtain 



(2) 



iedCr^'^Qp md^ (rjfsA'^^^A'f^^^Pi + vl^'^aAKaA^R) JaA ■ 



aA 



(24) 



Combining the contributions of Figs. ^ and |21 we can therefore write the total vertex 
function F'^ as 

F^ = ICda^Pqp [Vfs {^'(^aAKaA - KaA^'paA) 

aA 

(a1aA;;„A - ^lAKaA) Pr] JaA . (25) 

Clearly, if a = (3, i.e., we are talking about the diagonal vertex, the contribution vanishes, 
as is required by CPT invariance. 

3.2 CP-invariant limit 

Before proceeding further with the calculation, it would be instructive to check the form 
of the vertex in the CP-invariant limit. The CP transformation property of a fermion 
field ifj is given by 

(^^)7/;(x, t)(^^)-i = CCi^*i-x, t) , (26) 

where ^ can be called the CP phase of the field. For the scalar and pseudoscalar bilinears 
that occur in the Yukawa coupling of Eq. ()13|). the effect of CP transformation is then 
given by 

(^^)^1^2(^^)-' = 6?2*V^2jAl, 
Cr^)Vl75^2(^^)-' = -6^2*^275^1, (27) 

where the space-time points have not been explicitly mentioned. For all Dirac fermion 
fields, the CP phase can be absorbed into the definition of the antiparticle without any 
loss of generality. For Majorana neutrinos, however, that cannot be done because we have 
already defined through Eq. (fTTj) . In what follows, we will denote the CP phase of the 
field by ^a- 

We now take the CP transformation of the first term on the right side of Eq. (jlHj) . 

i^^)V^{vaKaAPL + Pr) rf„ C^^) = Cada{vAaAPR + Al^PL)/^aA^ . (28) 

If the Lagrangian is CP invariant, this CP transform should equal the hermitian conjugate 
term present in Eq. (|T3|) . This imposes the following conditions on the Yukawa couplings: 

A'* — ;^ n^A' 

KaA = CKaA- (29) 

If we put these conditions back into Eq. we find that the expression for the vertex 
function reduces to the form 

F^ = leaa^'qpVp Yl "^^^ \^PaAKaA (Pl - C/^Pr) + KaA^PaA {U^R " Pl)] JaA , (30) 

aA 
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where 

Ca/3 = CaVaCpV} = • (31) 

For Majorana fields, the quantity ^rj equals the CP eigenvalue of the 1-particle states 
[13, in Hi], so that denotes the relative CP eigenvalues of the two Majorana neutri- 
nos involved in the process. Thus, Eq. ()30|) shows that in the CP-conserving case, the 
transition amplitude is purely electric dipole type if the CP eigenvalues are the same for 
the two neutrinos, and purely magnetic dipole type if the CP eigenvalues are opposite. 
This is expected on quite general grounds 



3.3 CP-violating scenario 

We now go back to the general case, where the couplings do not necessarily obey Eq. 
and therefore the Lagrangian does not conserve CP. For such a case, we cannot use the 
vertex of Eq. ()30|) . Rather, we need to go back to the general form of the vertex given 
in Eq. ()25p. In fact, this expression can be written in a simpler form, noting the fact 
that the neutrino mixing matrix U and the phases rja are interconnected. In general, any 
mass matrix M can be diagonalized by a bi-unitary transformation involving two unitary 
matrices U and U': 

U'^MU = D , (32) 

where D is a diagonal matrix whose diagonal elements are the mass eigenvalues. Eq. (jH^ 
implies M^U'* = . For Majorana neutrinos, the mass matrix is symmetric. So, 
using M = and D = , we see that it is possible to choose U' = U*. Thus the 
bi-unitary transformation reduces to 

U^MU = D, (33) 

We will adopt this convention in what follows. This automatically implies |E| 

r]a = 1 for each a. (34) 

Of course, other choices are possible [T7j, but they all lead to same physical results. 
Therefore, the choice of Eq. ()34|) does not mean any loss of generality. 

To proceed, we make some simplifying assumptions. First, motivated by the fact that 
quark mixing is small, we will neglect it altogether and use the unit matrix in place of 

V. As mentioned earlier, the d-d'^ mixing is crucial for the terms that we have been 
considering in Figs. Hand|21 viz., terms having a factor of a down-type quark mass, let us 
assume that these are the only kind of mixings in the squark sector mass matrix. Further 
we assume that there is no intergenerational mixing in this sector, which means that the 

d squark of a certain generation mixes only with d^ of the same generation. Then the 
mixing matrix K in the squark sector is real, and is given by 





' c 


I) 




. -s 





(35) 



where in the final form, each entry denotes anxn diagonal matrix, where n is the number 
of generations, and 

C = diag (cos 6*1, cos 6*2, • ■ •)) 5" = diag (sin 6*1, sin 6*2, • • ■)• (36) 
More explicitly, we can write the elements of K and K as 

KkA = cos6k6k,A + sin 9k6k+n,A 

KkA = - sin 9k5k,A + cos 9k5k+n,A- (37) 

It should be recalled, from the definition of Eq. Q, that the index A runs over 2n values 
for n generations. 

The only intergenerational mixing in this picture is in the neutrino sector, which 
contains the phases of our concern. The vertex of Eq. ()25|1 can then be written as: 



(38) 



FPl + F*P 
where 

F = ^ TUd^ Bjk cos 9k sin 6k (^Jj,k+n - Jjk^ , (39) 

jk 

with Bjk defined by 

= X] (Ui'aUiP - UiaUi'p^X'ijkK'kj ■ (40) 
ii' 

Apparently, there is one big difference between this form of the vertex function and 
the expressions derived in an earlier work on the subject 7j- The expression in Ref. j7j 
contains two explicit factors of down quark masses. It should be emphasized that the 
extra factor, compared to our Eq. (j39|) . appears from an extra assumption. To appreciate 
this, consider what happens if all down-type squarks were really degenerate. The integrals 
defined in Eqs. p9|l and pn|) would have been independent of the squark index A in that 
case. Looking at Eq. we find that in Eq. the sum over A would have involved 
only the combination 

J2k:aK,a = 0, (41) 

A 

because of the orthogonality of two different rows of the unitary matrix K. The result 
implies that the vertex has to involve only differences of squark masses. These differences 
have to be small from various phenomenological requirements, and one often assumes 
that the differences are proportional to the corresponding quark masses. In particular, 
Eq. ()39|) shows that in absence of intergenerational mixing, the contribution from the two 
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squarks of a single generation vanishes if they are degenerate. Customarily, one assumes 
that the mass squared differences of the two squarks in the same generation is of the form 



Mm, 



(42) 



where ruq is the mass of the quark in the same generation and M is the average mass of 
the two squarks. It is this assumption which produced the extra factor of quark mass in 
the expression of Ref. [7j. We will carry on our discussion without this specific form for 
the squark mass differences. The phenomenological consequences of the decay rate, and 
the resulting bounds obtained thereby, have been discussed in earlier papers [ZllHl- In our 
discussion, we highlight CP violating features of this amplitude. 

4 Exploring phenomenological consequences 
4.1 The mass matrix 

Since the mixing matrix U is present in the vertex function through the combination 
of Eq. ()40|). the radiative decay amplitudes will clearly involve CP violating phases. To 
analyze how they appear in the amplitude, we first look into the generation of neutrino 
mass and mixing in this model. In the flavor basis for the neutrinos, the mass terms can 
be written in the form 



where the matrices appearing in the two chiral terms are related by the condition of 
hermiticity of their sum. Because of this relation, we can work entirely with just one 
part, say M, and call it the mass matrix. Diagonalization of this matrix was discussed in 
Eqs. ()32|) and (jH^ . Once this part is diagonalized and the eigenbasis is determined for 
left- and right-chiral fields, the other term is automatically diagonal in the same basis. 

It is to be noted that trilinear RPV couplings do not generate neutrino masses at 
the tree level. The masses are generated through 1-loop diagrams which look like the 
diagrams of Figs. Q and |21 without the photon line. If we calculate these mass diagrams 
with neutrino flavor states on the outer lines and with vanishing external momentum, we 
obtain 




(43) 




(44) 



with the integral J^^ defined as 




(45) 
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and the couplings defined by 

KaA = ^ Kjk^jaKlA , 
jk 

KaA = ^Kjk^kaKjA, (46) 
jk 

which are reminiscent of the couplings A and A in the neutrino flavor basis. The integral 
J' itself is logarithmically divergent, but we need to remember that the argument around 
Eq. (|4ip guarantees that Mjj/ vanishes if all squarks are degenerate. Since the divergence 
is independent of all masses, it cancels anyway. Of course, this has to happen since the 
theory is renormalizable and there is no tree-level mass term. 

In the simplified scenario advocated above where we neglect quark mixing and inter- 
generational squark mixing as well, the mass formula trivially reduces to 

M,i> = J2 {KjkK'kj + K'jkKkj) sin Ok cos 4 ['Jj,k+n - J'jk) > (47) 

jk 

which is what we will use below. 



4.2 Ansatz for a simplified analysis 



Phases can enter the mass matrix from any of the couplings A^^^. For three generations, 
there are 27 such couplings, and the phases are all independent. For illustrative purpose, 
we make a simplifying assumption that the phases can be written in the form 



A' 



ijk 



A' 



ijk 



(4J 



Since we are considering leptonic processes, it would have been most natural to focus 
on the phases Ui only, which pertain to lepton doublets. However, it is easy to see that if 
these are the only RPV couplings, i.e., the couplings are of the form 



A' 



ijk 



a: 



ijk 



(49) 



the model conserves CP. The reason is that in this case, we can absorb all phases of the 
X[ji.^s into the definition of the leptonic superfields Lj. The phases of these superfields will 
then reappear in the Yukawa couplings with the relevant Higgs doublet superfield, but 
those couplings can again be made real by readjusting the phases of the superfields Ej. 
Thus, at the end, no CP violating phase remains. If the phases denoted by (p and x 
non-zero, the possibility of CP violation remains, but the phases denoted by u continue 
to be irrelevant for the reason stated above. 

So we drop Ui altogether. Eq. PHjl now boils down to 



A' 



ijk 



A<t>j+Xk) 



\jk 



(50) 
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Through Eq. (jTfj) . this will introduce phases in the neutrino mass matrix. This mass 
matrix is diagonalized by the matrix U through Eq. ()33|). so that the phases in U are 
functions of the 0's and the x's- The functional relation is, however, non-trivial for three 
generations, where there are two Majorana-type and one Dirac-type phases in the neutrino 
mixing matrix |18j . 

We, however, note from the expressions in Eqs. (|39p and (j4Up that phases can enter 
the amplitude only through the quantities called Bj^. Clearly, the combinations occuring 
in Eq. (jH^ are very different from those occuring in the mass matrix, Eq. (^Tj), because 
the two equations involve different integrals and different combinations of the couplings. 
This implies that radiative decay amplitudes and neutrinoless double beta decay can in 
principle provide independent probes of CP violating phases. 

5 Conclusions 

Our primary intention has been the derivation of the radiative Majorana neutrino decay 
amplitude in the i?-parity violating supersymmetric framework, with a careful handling of 
the phases from different sources that enter the analysis. Our Eq. (^3)) is the most general 
expression for the radiative decay amplitude involving the A' couplings. In the derivation, 
we have not taken recourse to any specific assumption regarding the masses and mixings 
of fermions and sfermions. We have relied only on the general field-theoretic features, 
highlightling some subtle aspects regarding the Majorana character of the neutrinos which 
have not been clearly dealt with in the literature. We have shown, even in this general 
framework, that the radiative decay amplitudes depend on the quark masses and on the 
mass-squared differences of squarks. Our general formula can be used for all sorts of 
different realizations of fermion and sfermion masses and mixings. We have provided an 
illustrative example in Sec. 14.21 where we have pointed out that radiative decay amplitudes 
and neutrinoless double beta decay can in principle provide independent probes of CP 
violating phases. We do not discuss the numerical implications of our result, because they 
are of the same order of magnitude as those obtained in earlier papers, e.g., in Refs. [HIH]. 

As mentioned in the Introduction, we have not included the contributions from the 
A couplings. In fact, it is trivial to include these contributions. The diagrams will now 
involve charged leptons and sleptons in the internal lines, instead of quarks and squarks. 
Thus the contributions can be directly adapted from the A'-induced contributions by 
making suitable substitutions. 
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